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Abstract 

We discuss the properties of strange matter, in particular the minimum of 
the energy per baryon number as a function of the strangeness fraction. We 
utihze both the MIT bag model and the Color Dielectric Model and compare 
the energy per baryon with the masses of hyperons having the correspond- 
ing strangeness fraction, which are coherently calculated within both models. 
We also take into account the perturbative exchange of gluons. The results 
obtained in the two approaches allow to discuss the stability of strangelets. 
While the MIT bag model and the double minimum version of the Color Di- 
electric Model allow the existence of strangelets, the single minimum version 
of the Color Dielectric Model excludes this possibility. 



PACS: 24.85.+P, 25.75.-q, 14.20.Jn, 12.39.-x, 12.38.Mh 



1 Introduction 

It has been suggested long ago that hypermatter and/or strange quark matter 
can be realized in central relativistic heavy ion collisions, in which either multi- 
hypernuclear objects (strange hadronic matter) or strangelets (strange multiquark 
droplets) could be produced. The occurrence of the latter would be an unambiguous 
signature of the transient existence of a deconfined, strangeness rich, state of quark 
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gluon plasma (QGP). Of course the relevance of these objects highly depends on 
their stability after formation, and hence on the possibility of their detection. 

An enhancement of strange particle production has been observed since the early 
90's in many experiments with relativistic heavy ion collisions 0]. Indeed a large 
number of ss pairs can be produced in a single central event: the antiquarks s com- 
bine with the abundantly available light quarks u and d, thus producing antikaons, 
which rapidly leave the fireball region. The residual system turns out then to be 
a strangeness rich matter which hadronizes with a copious production of strange 
particles, especially K mesons and A hyperons: strongly enhanced yields have been 
by now assessed [Q . Yet this outcome cannot be considered as a reliable signature of 
quark-gluon plasma (QGP) formation, since kaons and hyperons can be produced 
in hadronic reactions as well, before the nuclear fireball reaches equilibrium 0. 

The case would be different if, after formation of the deconfined plasma, this 
strangeness rich matter could coalesce into colorless multiquark states, the so-called 
strangelets. The prompt anti-kaon (and also pion) emission from the surface of the 
fireball could, in addition, rapidly cool the QGP, thus favouring the condensation 
into metastable or stable droplets of strange quark matter 

Long after the first definition of strangelets, by Bodmer and Chin and Ker- 
man |^, Witten conjectured [§] that strange quark matter can be absolutely stable 
(namely more stable than ordinary nuclei), a result obtained within the MIT bag 
model with parameter 5 ~ 58 MeV fm~^. The main reason underlying this situ- 
ation is the lowering of the Fermi energy introduced by the new degree of freedom 
(the s-quark), which in turn lowers the global minimum of the energy per baryon 
{E/Ab) with respect to ordinary nuclear matter. 

It is worth mentioning that strangelets can also be produced by coalescence of 
hyperons ||^: it consists of the fusion of a few hyperons and nucleons and does not 
imply the existence of the QGP phase: however this process can typically produce 
only small baryon numbers (e.g. the H dibaryon). Hence the detection of heavy 
strangelets {Ab > 20 — 30) should remain an unambiguous signature of QGP for- 
mation. For a review on the properties of strange matter and strangelets we refer 



the reader to the papers by Greiner and Schaffner [Ty, [TT|. 

Up to now the properties and stability of strangelets have been discussed within 
the MIT bag model or, similarly, a Fermi gas model stabilized by the vacuum pres- 
sure B] the pioneering work by Fahri and Jaffe []T2| also includes 0{as) corrections 
to the properties of bulk strange matter and finds the stability conditions: according 
to these authors heavy, slightly positively charged, strangelets could be more stable 
than ordinary nuclei. Greiner et al. find (in pure MIT bag, without color ex- 
change contributions) that also light strangelets (with Ab < 20) are metastable due 
to finite size and shell effects. A detailed calculation of strangelet properties within 
the MIT bag model, including shell effects and all the hadronic decay channels has 
been performed by J. Schaffner et al. [0|: a valley of stability clearly appears for 
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Ab = 5 16 with charge fraction Z/A between and —0.5. On the other hand, 
strangelets having a larger mass should be positively charged according to the results 



of Ref. 



In this paper we want to discuss the properties of strange matter, in particular the 
minimum of the energy per baryon number as a function of the strangeness fraction. 
We utilize both the MIT bag model and the Color Dielectric Model (CDM) and 
compare the equilibrium energy of the strange matter with the masses of hyperons 
having the same strangeness fraction; the latter are coherently calculated within 
both the MIT and CDM models. The main goal is to find out whether and to which 
extent the stability of strange matter and/or strangelets with respect to ordinary 
hadronic matter depends on the model employed to describe the confined system of 
quarks. 

We consider homogeneous quark matter made up of u, d and s quarks. We 
assume that, during a high energy collision between heavy ions, this state of matter, 
if formed, can only survive for a very short time, so that it has no time to reach 
jS equilibrium; hence we do not impose chemical equilibrium on the density of the 
strange quarks, limiting ourselves to assume that there exists in the system a strange 
fraction Rg = Ps/p, p being the total baryon density of quarks and ps the baryonic 
density of strange quarks. Our aim is to find out whether this system is more or 
less stable than hyperons, in order to understand which state is more likely to be 
produced in heavy ion collisions, either hyperons and strange mesons or strangelets. 

The energy per baryon number in the mean field approximation is a function of 
the baryon density, p, and of Rs- Thus we consider the surface which represents 
the energy per particle, E /Ab, versus p and Rg and concentrate on its sections with 
constant Rs'. the resulting curves represent the energy per baryon number, at a 
fixed Rsi as a function of p, and in general they present a minimum. The minimal 
energy per baryon number can then be studied as a function of the strange fraction 
Rg. In order to describe the system of three flavors we shall employ two different 
models: the MIT bag model and the CDM. We also consider the effect induced 
by the introduction of perturbative gluons. Electromagnetic interaction has been 
neglected in this paper and therefore the minimum of the energy corresponds to 
an equal number of u and d quarks. In the MIT bag model, u and d quarks are 
massless, while in the CDM they are massive: thus it will be interesting to compare 
the results obtained in each case, in order to find out the role of light quark masses 
on the bulk properties of strange quark matter. Another difference between the two 
models is that in the MIT bag model there is a sharp transition from the inner region 
to the outer region of the nucleon. In the CDM, on the contrary, both a two phase 
and a single phase scenario can be obtained. As we will see, these two possibilities 
give rather different results for the stability of strangelets. 

Since we consider an infinite and homogeneous system, while strangelets are 
finite objects, one should be careful in drawing conclusions about the stability of 
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strangelets on the basis of strange matter stability; however the energy of the infinite 
system appears to be a lower limit with respect to the envelop of strangelet energies 



versus strangeness fraction: the latter was nicely illustrated by Schaffner et al. 
calculating the strangelet masses within the MIT bag model with shell mode filling. 
We report their result in Fig. |l], which will be useful for further comparison with 
our results. Here we simply recall that surface effects, which we do not consider, 
would increase the energy curves of bulk matter, typically of 50-100 MeV: hence, if 
hyperons should turn out to be more stable than strange matter, then this would 
exclude also the stability of strangelets. If, on the contrary, strange matter is more 
stable, then this provides only an indication in favour of stable strangelets. 



2 Strangelets in the MIT bag model 

The MIT bag model has been widely utilized in the past, both for strange matter 
and for strangelets. For detailed derivations we refer the reader to the literature 



12, 13, 16 



2.1 MIT bag model without gluons 

We shall use here the simplest version of the MIT bag model, not including one 
gluon exchange corrections. Therefore, we have two model parameters: the vacuum 
pressure B and the strange quark mass m^. We have used different values of these 
parameters, in order to discuss various possible scenarios. 

Calculations of hadron spectroscopy [|l^ indicate for B a value of the order of 60 



MeV/fm^. However this value is generally not adopted in quark matter calculations, 
because it produces too much binding (notice that it is close to the 58 MeV/fm^ of 
Witten's conjecture); calculations of the hadronic structure functions suggest 
a value of -B ~ 100 MeV/fm^. Hence, besides -B = 60 MeV/fm^, we consider also 
B = 100 MeV/fm^ and an additional value, B = 150 MeV/fm^. The latter does not 
correspond to calculations of physical quantities, but it has been indicated in the 
literature |T^ as a sensible value in a comparison between the results of MIT bag 
model and lattice QCD at finite temperature. 

The single flavor contribution to the energy density of the system is given by: 



/dk 
-—,Ef (k) 9 {kp^. - k) , 



(2vr) 

where Ef{k) = ^Jk"^ + mj- and kpj, is the Fermi momentum of flavor /. We use 
natural units {h = c = 1). Since u and d quarks are massless, their contribution to 
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the energy density is: 



(2vr) 



(2) 



while the strange quark contribution reads: 



mt In 



+ ki 



kps + \ k% + ml 



kl^ + m2 (24^ + 



(3) 



being the mass of the quark s. 

The total energy density of our system turns then out to be: 

^tot = '^(-uA + es + B . 



(4) 



The dependence of the above formula on Rg and p can be easily found by recalling 
that: 

Ps = Rsp (5) 

and: 



kf 



kp 



2„ nV3 



37r2 



Pil-Rs 



1/3 



(6) 
(7) 



p being the total baryon number density in the system (p = Ab/V). In the above 
the color degeneracy and baryon number 1/3 of the quarks have been taken into 
account. 

From the above formulas we calculate the energy per baryon number to be: 



tot 



^tot 

p ' 



In Fig. 1^ the results of the minimal energy per baryon (^) corresponding to -B = 
60, 100, 150 MeV/fm^ are shown as a function of Rg. For each value of B we explore 
three different values of the strange mass, rris = 100, 200, 300 MeV, and we compare 
these results with the experimental nucleon and hyperon masses (full circles). We 
have also evaluated, according to formula (3.6) of Ref. JlBl, the baryonic masses 
which are obtained within the same model employed for bulk strange matter, using 
the same sets of bag parameter and strange quark mass. Here we set to zero the 
perturbative one gluon exchange corrections, which are, instead, taken into account 



by De Grand et al. [0] and which will be considered in the next subsection. For 
baryons, the fraction Rg is assumed to correspond to the fraction of s quarks with 
respect to u and d quarks in the considered hadron. 
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As it appears from the figure, the three hnes corresponding to the different values 
of TUs are much lower than the experimental hyperon masses for B = 60 MeV/fm^ 
and B = 100 MeV/fm^, while this is not the case for B = 150 MeV/fm^ and 
ms=300 MeV. In this instance the experimental hyperon masses are lower than the 
bulk matter energy. However, if we compare the energy of strange matter with 
the corresponding theoretical masses of the various hyperons, we find that strange 
matter is always lower in energy, and thus more stable. In agreement with the results 
of Ref. [0, which are shown in Fig. ||, we can conclude that the MIT bag model 



without perturbative gluon corrections allows the existence of strangelets. It might 
be interesting to notice also that in the MIT bag model a minimum at a finite value 
of Rg is always present. The only parameter set not showing a minimum corresponds 
to 5 = 60 MeV/fm^ and m, = 300 MeV. 



2.2 MIT bag model with perturbative gluons. 

We consider now the effects of the introduction of perturbative gluons in the cal- 
culation. At first order in a^, two contributions to the energy can be considered, 
the direct and the exchange one. Since the system is globally colorless the direct 
term vanishes, while the exchange one gives the following contribution to the energy 
density of quarks of flavor / |l^ :Q 



^OGE 



Here: 



as 4 I 4 3 
-nif < Xf 



3i 2 1 
+- In^ — 

2 \Vf 



In 



3 In 



Vf 



m 



fVf 



^fVf 



+ 



[r]fXf - In {xf + r]f)] 



(9) 



Vf 



rrii 



Xt 



and /i is a renormalization scale, for which we choose the value /z = 313 MeV, 
according to Ref. In Fig. ^ we show e^'^^ as a function of hp for various values 
of m. It is evident that, for small values of m, the contribution is always repulsive, 
while for m > 200 MeV there exists a range of densities for which it is attractive. 
Hence the effect of including perturbative gluons in the energy density of the system 
will crucially depend on the fraction of strange quarks, as well as on the value of 

^ Actually the third term in square brackets has a different sign, with respect to Ref. in 
agreement with a recent re-derivation of the formula |l9[. 
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mg. If the quark mass vanishes (mj = 0) eq. (j^) reduces to: 



^^" = 54.,,. (10) 

Yet in the following we shall use a small nonzero mass for the u and d quarks, 
^u,d = 4 MeV. For sake of illustration, we adopted two different values for a^, a 
small perturbative value (a^ = 0.5) and the canonical value which was employed by 
DeGrand et al. (a^ = 2.2), to reproduce the hyperon masses. The corresponding 
results are illustrated in Figs. § and |. 

By comparing Fig. |^ and ^ we can see that the gluon effect is always repulsive 
at low while at large Rg it is repulsive for = 100, 200 MeV, and attractive 
for iris = 300 MeV. This is even more evident in Fig. ^, where the effect of gluons is 
stronger, due to the larger value of a^: in this case the attractive effect at large Rs 
is so important that the curve corresponding to mg = 300 MeV becomes the lowest 
one, in contrast with the situation shown in Figs. |^ and |[ 

From Fig. § we can see that, even after the inclusion of perturbative gluons, 
strangelets are more stable than hyperons for almost all values of the model param- 
eters. However, when we use the stronger coupling of Fig. ^ the stability of strange 
matter (and hence strangelets) becomes questionable, particularly for low values of 
the strange mass m^. Only for = 300 MeV the theoretical masses of hyperons 
always lie above the energy of bulk matter (not so the experimental masses). It 
is worth keeping in mind that among all the situations illustrated here, only the 
upper left panel of Fig. |^ utilizes parameters close to the ones of DeGrand et al. |T^ 
{B = 59 MeV/fm^, m, = 279 MeV, = 2.2): indeed the stars (m^ = 300 MeV) re- 
produce fairly well the experimental masses, and for this choice of parameter values 
the bulk strange matter turns out to be favoured. 

From this analysis we can conclude (in agreement with previous findings) that, 
apart from rather extreme choices of the model parameters, metastable strangelets 
can exist in the MIT bag model. 



3 Strangelets in the Color Dielectric Model 



The Color Dielectric Model provides absolute confinement of quarks through their 
interaction with a scalar field x which represents a multi-gluon state and produces 
a density dependent constituent mass (see for example the review articles [^ - [p^ ). 



Several versions of CDM have been employed to calculate hadronic properties 
| 29[| and quark matter equation of state [|30[ 



m 



The typical Lagrangian of the CDM is the following [31 



f=u,d,s 



X 



f=u,d 
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+ 2 id,X? -U{x)- (X) F;^F'^^^ , (11) 

where V'/ are the quark fields, A° is the (effective) gluon field, F^^ its strenght tensor 
and X is the color dielectric field. 

The u and d quark masses are the result of their interaction with the x^field and 
read: 

= — (12) 
X 

where is a parameter of the model and /^r the pion decay constant, which is fixed 
to its experimental value, /^r = 93 MeV.| For the strange quark mass we consider 
two different versions of the 3-flavors CDM, namely a scaling model, with 

gU + Am g'f^ 

ms = = , (13) 

X X 

and a non-scaling model, with a constant shift of the s-mass with respect to the 
u, d-one: 

= + Am = m,u d + Am . (14) 
X 

In the above g' and/or Am is another parameter of the model. 

Concerning the color dielectric field, there exist several versions, both for its 
coupling to the gluon tensor and for the potential U{x)- We adopt here both the 
single minimum (SM), quadratic potential: 

f/sM(x) = ^MV, (15) 

which introduces the third parameter of the model, M (the mass of the glueball), 
and the double minimum (DM), quartic potential: 

The latter introduces an extra parameter, the bag pressure B, while the parameter 
Xo is used to make the ratio xIxq dimensionless. The color-dielectric function, 
ft(x), is usually assumed to be a quadratic or quartic function of %: we will use both 
options and hence we set: 

X-' 



k{x) = [ — ] , with /3 = 2,4. (17) 



^ In some chiral invariant versions of the CDM the mass term is also coupled to the usual (ct, tt ) 
fields |||9[ ||. 
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The field equations will be solved in the mean field approximation and with- 
out gluons: the latter are subsequently taken into account as a perturbation. For 
homogeneous quark matter the color dielectric field must obey the equation: 



dU{x) 



dx 



x=x 



f=u,d 



)fn, T I . 

X 



where the second term on the r.h.s. will contain g or g', depending on the choice 
of the non-scaling model {g) or of the scaling one {g'), respectively. In the above 
equation: 



(ijfijf) = pi ix) = 6 



dk 



^fix) 



(27r) 



-.6 {kpf - k) 



(19) 



is the scalar density of quarks /. 

The unperturbed (i.e. without gluon contribution) energy density reads: 



— - < mt In 
87r2 I f 



f=u,d,s 



kp, + ^|k\ +m) 



+kFf yjkl^+m) {2kl^ + mj) } + f/ (x) 



(20) 



the quark masses being given by eqs. ([T2|) and (|r^) [or (|T4|)] with x = X- 

Beyond eo we have perturbatively taken into account, to order a^, the exchange of 
gluons, whose contribution to the energy density of an infinite, color singlet system 
is, in analogy with eq. (pi): 



(-OGE 



OGE 



V 



E 6 

f=u,d,s 



dq 



■VoGE{q)0 {kF, - k) 



(21) 



f=u,d,s 



4 J 4 



+ _3in 
2 \7]fJ \mf7]f 



2 



In 



Vf 



^fVf 



[Xfr]f - In {xf + r]f)] 



(22) 



Here the notations are the same as in eq. (^, but for the CDM definition of the 
constituent quark masses (in MFA); the effective strong coupling constant, according 
to the CDM Lagrangian (0), reads: 



Xo 
X 



(23) 
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and becomes very large at small densities, where x 0. As already remarked 
above, eq. (|21|) only contains the exchange term of OGE, the direct one vanishing for 
infinite quark matter: at small baryonic densities the attractive electric contribution 
dominates the energy density, while the repulsive magnetic contribution becomes the 
dominant one at large densities. 

Indeed the divergent behavior of the electric term for p ^ could prevent the 
reliability of a perturbative treatment of OGE. We have overcome this difficulty by 
taking into account the Debye screening of the gluon propagator in the presence of 
a polarized medium. This can be achieved by replacing (^) with the new effective 
coupling: 



a 



eff 



a. 



g{y) being the static limit of the polarization propagator 134 



1 



(24) 



(25) 



Actually this expression should be utilized in Voge{(1) of eq. ( pTj ) and again 
integrated. For simplicity, since the g-integration is extended only up to kpj, and 
the function g{y) varies at most by 9% in the range < y < 1, we have adopted 
q = kpj. 

The flavor dependent effective coupling reads therefore: 



a 



eff 



Xo 
X 



k\ 



iXo/x) rnjkl g{ 



(26) 



We notice that, taking into account the dependence of x upon the density, this new 



effective couphng vanishes at small densities: 



a 



eff 



while it goes to zero as k^ 



at high densities. At variance with the Debye screening in electrodynamics, which 
is mainly relevant at large densities, in the CDM the screening is large at small 
densities too, due to the divergence of 5,, which enhances the effect of medium 
polarization even at small densities. With these ingredients we shall compare our 
results for the minimal energy per baryon number with the hyperon masses, as 
they have been evaluated in the CDM. There exist in the literature two distinct 
calculations of this type and we shall consider both cases. 



3.1 Stability of strangelets in the CDM: I 

We consider here the work by Aoki et al. ||2^: these authors solve self-consistently 
the mean field equations for quarks, color dielectric field and gluons, starting from 
a CDM Lagrangian with the Double Minimum potential (p!6| ) for the color dielectric 
field. Concerning the color dielectric function, Aoki et al. choose /3 = 2 in eq. (0). 
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In order to reproduce the masses of the octet and decuplet baryons (in particular 
hyperons) the authors of Ref. ||26| employ both the scaling model and the non-scaling 
one. In addition, two different sets of parameters are used, whose values are dictated 
by two different and extreme choices for the "bag" parameter B: B^^^ = MeV, 
with two degenerate vacua, and a large bag pressure, B^^^ = 103.5 MeV. The latter 
value of B is chosen to be as large as possible, but with the requirement that the 
two-phase picture must hold inside hadrons. In their calculation only the strange 
quark mass has to be considered as a truly free parameter, the remaining ones having 
been fixed in a previous work on the non-strange baryons pSf. We only perform 



calculations with the second set of parameters, corresponding to the non-zero value 
of the bag pressure, since B = does not give quark confinement at low densities. 
The field equation for the color dielectric field, eq. (|TB|), now becomes: 



for the non-scaling model, and: 

(28) 

for the scaling model, respectively. The parameters of Ref . (with i? 4 = 103.5 MeV) 
appropriately converted to the notation of the present work, are the following 
(«, = (72/47r): 



c/ = 43.7MeV M = 1177MeV xo = 47.1MeV a, = 2.0 
Hence, we consider cases B and D of the work of Aoki et ai: 
B: scaling model, g' = 138.9 MeV; 
D: non-scaling model. Am = 212 MeV. 

We report in Table I the masses of a few hyperons (the ones of interest for the 
present work) as obtained in Ref. |2^, together with their experimental values: the 
agreement is remarkable. 

We evaluate the minimum energy per baryon number using cases B and D, 
both without and with the perturbative correction of the effective gluon exchange.0 

■^In the Appendix we discuss the existence of muhiple solutions for field equations when the 
potential U (x) has a double minimum. 
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Hyperon 


N 


A 




Q 







0.33 


0.67 


1 


rriexp (MeV) 


938 


1115 


1314 


1672 


niB (MeV) 


938 


1161 


1346 


1639 


niD (MeV) 


938 


1113 


1307 


1671 



Table I. Strangeness fraction, experimental masses and theoretical masses of hyperons, according 

to the calculation of Ref. . 



In Fig. g we compare our curves with the masses calculated in Ref. [^], which 



always include, as already mentioned, the gluon field correction in a self-consistent 
approach. In addition to the theoretical masses of hyperons, we also show in the 
figure their experimental values (represented by full circles) which, according to 
Table I, are always rather close to the calculated ones. 

As we can see from Fig. |^, the effect of perturbative gluons in this model is very 
small, due to the Debye screening, which we have included. Whether or not we take 
into account gluon corrections, strange matter always appears to be more stable 
than baryons. 

Hence from the cases considered here we can conclude that this version of the 
CDM seems to favour strangelets as a metastable form of matter. This is due to 
the fact that when the DM potential is used to study hadrons, i.e. confined objects, 
a large contribution to the hadronic mass is given by the space fiuctuations of the 
fields. When this version of the model is used to describe infinite quark matter, 
these contributions vanish due to the homogeneity of the system. For this reason, 
deconfined matter is favoured in this version of the model. This finding is consistent 
with the results of Refs. |]30| , It is also worth mentioning that, as shown in 



Fig. ^ the minimum of the energy per baryon number corresponds to i?^ = 0, at 
variance with the results of the MIT bag model. 

3.2 Stability of strangelets in the CDM: II 



In this section we follow the approach of J. McGovern |2^. The model Lagrangian 
is still the one reported in eq. (|rT|) with /? = 4 in the color dielectric function. 
McGovern uses only the scaling model and the Single Minimum potential, with 
different values of the parameters. In this case the behavior of a';/^{x) is even more 
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Hyper on 


N 


A 




Q 


Rs 





0.33 


0.67 


1 


rriexp (MeV) 


938 


1115 


1314 


1672 


m,et I (MeV) 


1227 


1427 


1638 


1866 


mset II (MeV) 


1268 


1417 


1691 


1928 



Table II. Strangeness fraction, experimental masses and theoretical masses of hyperons, as 

calculated in Ref. p7| . 



divergent, for small densities, than in the case /3 = 2 previously discussed: hence 
the use of Debye screening in the effective strong coupling constant is mandatory. 
In Ref. [0 two different sets for the model parameters are used: they allow 



to satisfactorily reproduce the splittings between hyperon masses, but the absolute 
values of the masses themselves are generally too large. The latter are reported in 
Table II. 

In this model, the field equation for the scalar field becomes: 

M'x = # [P's ix) + Pi m + ^-^P's ix) . (29) 
X X 



The parameter sets used in Ref. are: 
Parameter set I: 

g = 46.7 MeV M = 2354 MeV xo = 20 MeV a, = 0.3380 



Parameter set II: 

g = 107.527 MeV M = 1000 MeV Xo = 45.4 MeV a, = 0.3533 



We fix the ratio: g'/g = 1.58, in agreement with Ref. ||27 



As in the previous section, in Fig. |^ we compare our results for the minimum en- 
ergy per baryon number both with the experimental and the theoretical masses. As 
we can notice, also in this case the inclusion of perturbative gluons practically does 
not alter the quark matter energy, since here the Debye screening is overwhelming. 
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Strange matter turns out to be well above the experimental hyperon masses and 
below the theoretical ones. We must remember that in our calculations we don't 
consider surface effects, which would increase the energy density of a finite system 
of 50 -7- 100 MeV, as we can deduce by comparing our results in the MIT bag model 
with Fig. m. In this perspective, only for i?^ ~ | strangelets are (marginally) allowed 
by the present calculation and a more refined one, taking into account surface energy 
contributions, could clarify the situation. 

We can therefore conclude that, but for Rg ~ |, the existence of metastable 
strangelets is excluded in the Single Minimum version of the CDM. 



4 Conclusions 

In this work we have computed the energy per baryon number of infinite quark 
matter having a fixed strangeness content. We have compared this quantity to the 
mass of hyperons having the same strangeness fraction and calculated within the 
same model and for the same parameter values that we adopt in our calculations. 

Our analysis shows that the existence of strangelets is supported only by those 
models which entail a two-phase picture of hadrons, namely which mantain a false 
vacuum inside hadrons. This happens both in the MIT bag model, and in the Double 
Minimum version of the Color Dielectric Model: as we have seen in Sections 2 and 
3.1, in these cases the minimum energy per baryon number versus the strangeness 
fraction Rg turns out to be lower than the masses of hyperons with the same Rg. 

As we have seen in Figs. H, ^, ^, the mass gap between hyperons and strangelets 
can be as large as 300 MeV, and the metastability of strangelets would be confirmed 
even by taking into account surface energy effects, typically of the order of 100 MeV, 



in agreement with the findings of Ref. 



On the contrary, the Single Minimum version of the CDM does not support the 
existence of strangelets: independently of the set of parameter values, the computed 
masses of hyperons are larger than the corresponding strange matter energy by 
about 50 ^ 100 MeV. Only for Rs = ^ the gap is 150 MeV and we cannot exclude 
completely the presence of strangelets with this strangeness fraction. 

This outcome points out that the stability of strangelets appears to depend rather 
crucially on the model employed: CDM supports stable strangelets only in the DM 
version, but not in the SM one. This fact could set serious challenges to the search 
for strangelets in relativistic heavy ion collisions, as a signature for the quark gluon 
plasma phase, out of which strangelets could be formed. 
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5 Appendix 



There can be multiple solutions of field equations, when for the dielectric field x 
a potential having two minima is used. This possibility is known since the first 
calculations within the color-dielectric model ||35|, ^ . 

In the above quoted references, the solutions of the field equations were non- 
topological solitons. In the present work we are discussing plane-wave solutions for 
the quarks, which are considerably simpler. The solutions of the field equations for x 
are graphically shown in Fig. |], where they appear as the intersection of the dotted 
line with the dashed (at lower density) or continuous lines. As it appears from 
Fig. 1^, when the density is low enough there can be three solutions. The solutions 
are characterized by their behavior at low density. One solution corresponds to 
X — >■ as p — > 0. The other two solutions correspond to x close to the value of 
relative maximum of the potential or to the relative minimum, respectively. The 
solution in which x is near the relative maximum is unstable, and it corresponds to 



the unstable solitonic solution pointed out in Ref. |]36[ . 

At low density there are therefore two stable solutions. The "true" solution is 
clearly the one of minimal energy. As it can be seen in Fig. ^, but for very small 
value of the density, the solution of minimal energy is the one in which x is close to 
the value of the relative minimum of the potential. 
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Figure 1: The energy per baryon of all possible strangelets with As < 40, for a bag 
constant B^^'^ = 170 MeV versus the strangeness fraction fg = 3Rs = \S\/Ab- The 
solid line connects the masses of the nucleon and of the first lowest hyperons: it 
represents free baryonic matter. (Taken from Ref. []I3| ) 
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Figure 2: Minimal energy per baryon number in the MIT bag model, as a function 
of the strangeness fraction Rg = Ps/p, for various choices of the values of the model 
parameters. The continuous line corresponds to = 100 MeV, the dashed hne to 
nis — 200 MeV and the dotted hne to = 300 MeV. Pull circles correspond to 
experimental masses, the other points to the masses evaluated in the model, with 
THs = 100 MeV (open triangles), = 200 MeV (fuU triangles), = 300 MeV 
(stars), respectively. 
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Figure 3: OGE contribution to the energy density divided by q;^, e'^'^^/ct^, as a 
function oi kp for different values of the quark mass. 
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Figure 4: Minimal energy per baryon number in the MIT bag model, including the 
OGE potential with = 0.5, as a function of the strangeness fraction Rg = Ps/p- 
The continuous line corresponds to = 100 MeV, the dashed line to rUs = 200 MeV 
and the dotted line to rris = 300 MeV. Full circles represent the experimental masses, 
the other points refer to the masses evaluated in the model, with = 100 MeV 
(open triangles), = 200 MeV (full triangles), = 300 MeV (stars), respectively. 
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Figure 5: The same as in Fig. 0, but for as = 2.2. 
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Figure 6: Minimal energy per baryon number in the CDM, as a function of Rg = 
Ps/ p, for the cases B and D with and without gluons. Full circles are the experimental 
hyperon masses, while triangular dots are the masses calculated in Ref. p6| . 
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Figure 7: Minimal energy per baryon number as a function of Rg = Ps/ p for the Sin- 
gle Minimum version of the CDM. The various panels correspond to: (a) parameter 
set I without gluons, (b) parameter set I with gluons, (c) parameter set II without 
gluons and (d) parameter set II with gluons. Full circles are the experimental baryon 
masses, while triangular dots are the masses calculated in Ref. 
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Figure 8: Solutions of the field equation for the scalar field, for two different values 
of p: the dotted line corresponds to df/^Af/dx, the dashed and continuous lines 
correspond to the r.h.s.of eq. (^), assuming p = 0.01 fm~^ (dashed line) and p = 
0.05 fm~^ (continuous line), with fixed Rg = 0.05 in both cases. 
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Figure 9: Total energy per baryon number as a function of p for Rg — 0.5: the solid 
line corresponds to the curve calculated with x near to the value of the relative 
minimum of the potential, the dashed line corresponds to x near to the absolute 
minimum of the potential. 
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